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We have studied the link-integration method for the improved actions. With this method the r\ parameter in 
the medium to strong coupling regions is obtained. Effects of the self-energy terms for the rj parameters are small 
in the regions of (3 and rj studied. After these investigations, the anisotropic lattice is used for the calculation of 
transport coefficients of the quark gluon plasma. 



1. Link Integration For Improved Action 

If R is an external source field for link variable 
U, the expectation value of U is given by, 
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Z(R) is expressed by the modified Bessel function 
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Q = 2Re(det(R)) 

P{x) = l + xTr{RRt) 

+ \x 2 [(Tr(RR^)) 2 - Tr((RR^) 2 )] 
+x 3 det(RR^) 
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Similarly dZ(R)/dR is written by the I\ and 

The path of the integration is a closed circle on 
the complex plane x. In principle its radius r is 
arbitrary, but numerical integration requires ade- 
quate radius. We apply Simpson method for the 
numerical integration, and search for the region 
of r and number of the division N where < U > is 
stable with the change of r and N. It is observed 
that the < U > depends strongly on r and N. 
For example if we take N = 50, there appear spu- 
rious plateaus, which disappears with increasing 
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Figure 1. An example of r%* x , , min 
a fully thermalized configuration 



N . However there is a region of r in which < U > 
is stable for the change of r and N. It becomes a 
little wider as N is increased^]. This region of r is 
called an optimal region of integration;r opt . The 
link integration should be done with r opt 

The r opt are shown as a function of (3 in Fig.0. 
Notice that r opt region changes with the position 
of link variables on a configuration and also with 
configurations due to the fluctuation of gauge 



1 These features are seen not only in the asymptotic ex- 
pansion of the modified Bessel functions but also if Taylor 
expansions of them are applied. 
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fields. The results shown in Fig|l] are not the av- 
erage over the fluctuations. But the fluctuation 
of the r opt region is not large. If we choose, 
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it has been in an optimal region of r throughout 
the link variables and configurations. The r opt for 
space link U space is larger than that of Uu m e when 
£ > 1; r° p l ce > r°[ me , and they become smaller 
for larger j3. 

In the case of improved actions, the number 
of link U which are simultaneously integrated 
in a loop becomes much smaller than the case 
of standard action. Therefore the link integra- 
tion method is not effective for the calculation 
of smaller Wilson loops. The suppression of the 
fluctuation is impressive for W(6 x 6) but not 
for W (4 x 4) for /3 — 2.5. However for the cal- 
culation of T] in the smaller (3 region, the use of 
the link integration method has been indispens- 
able. We have used it for the calculation of r\ at 
(3 = 3.5, 3.05, 2.5 for Iwasaki action. 

2. Self-energy Effects on r\ 

Lattice potential is defined by the ratio of Wil- 
W(p r) 

son loops, V(r) = log(— — — — -). The (, R is 
W(p + 1, r) 

determined by the matching of the potential in 
space and temporal direction. 



V s (Z B ,r) = Vt(£ B ,t = Z R xr) 
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The 77 parameter is defined by r\ = £r/£b How- 
ever, V s and V t includes self-energy contributions 
which may be written as, 

V a fo ) r) = V?fo) + V„fo,r) (6) 

Similarly for V t (£ B ,t). 
Due to the anisotropy, V°(£ B ) ^ V t °(£ B ). 
For the standard action, the effect of the self- 
energy term on 77 has been studied by Bielefeld 
groupQ. It is reported that the effect on 77 is 
~ 0.01. In this report we study its effects for 
Iwasaki's improved action S. In order to get rid 
of the effect of self-energy term V° , we employ a 
subtraction method, 



Vf ub tt B ,r) = V s (H B ,r)-V s (Z B ,r a ) 



(7) 
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Figure 2. 77 s " 6 (filled circle) and v Nosub without 
subtraction(open Symbols) at £r = 2.0 



similarly for V t Sub . The to and ro should satisfy 
to = £,n r o and we take ro 
matching for y s " h [|§. 



3. Then we apply the 



Vf vb fe ) r) = V t ™fe t t = Z R xr) 



Sub, 
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We denote the 77 determined in this way as r) Sub . 
In Fig.0 we have shown our results. It is found 



that the differences of the 77 



Nosub 



and 77 



Sub 



are 



small; < 0.015. 

Combined with the results of Bielefeld group 
for the standard action, we conclude that the cor- 
rection to the global dependences of the 77 param- 
eters for the class of improved actions(Symanzik, 
Iwasaki etc) which we have presented at lattice 
99 § is really small 

3. Quark Gluon Plasma on Anisotropic 
Lattice 

As applications of the anisotropic lattice, we 
have started a simulation of transport coeffi- 
cients of the quark gluon plasma and a study of 
the heavy quark spectroscopy 0). 

We reported the transport coefficients from lat- 
tice simulations on isotropic lattice||||. Our 
results have been impressive and encouraging 
for the phcnomcnological study of the quark 
gluon plasma, in the sense that they are located 
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close to the simple extrapolation of the pertur- 
bative results on the figure at high tempera- 
ture limitj^]. However they depend on the ansatz 
for the spectral function of the Matsubara green 
function(G Atl y) of energy momentum tensor. The 
aim of this work is to improve the resolution of the 
G^ v by using the anisotropic lattice, and check 
the ansatz of the spectral function. 

The calculation of the transport coefficients of 
quark gluon plasma are formulated in the linear 
response theory. They are calculated by the G^ v 
of energy momentum tensor. For the pure gauge 
models it is written as, T^ v — 2 Tr [F^F^ — 
jfifivF p(J F pcr } On an anisotropic lattice, the field 
strength tensors are written as follows, 



G„. p=3.3 



alF? j FZ = p-ri(3-TrU i j) 



a a F i4 F i4 



TrU, 



(9) 
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We have started from test run of calculating G pv 
on 8 3 x 8 lattice with Iwasaki action at (3 = 3.3 
and = 2.0. As a check of our calculations on 
the anisotropic lattice, we have compared it with 
the one from an isotropic lattice 8 3 x 4, for which 
the scale in temperature direction is changed to 
in = 2 lattice. They are shown in Fig.||. 

It is seen that the both green functions co- 
incide with each other. The merit of applying 
an anisotropic lattice is that we could determine 
the transport coefficients under the ansatz of the 
spectral function on small lattices. 
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Because 3 independent parameters are deter- 
mined by the 3 independent data points of G\i 
on an anisotropic lattice. 

We have proceeded to the calculation of Mat- 
subara green function on 24 3 x 16 lattice with 
£r = 2.0. At this point, the number of data 
is not large enough to determine the G pu . But 
we expect that the check of the ansatz given by 
Eq.[n], could be done on the anisotropic lattice. 
In addition we are planing to determine spec- 
tral functions by the maximum entropy method. 

2 Notice that perturbative formula breaks down around T c 
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Figure 3. Matsubara green function G12 on 
8 3 x 8 lattice with £g = 2.0(open circle), and on 
isotropic 8 3 x 4 lattice(ffiled square) at (3 
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They will be reported in the forthcoming publi- 
cations 
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